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Abstract 

Using similarity transformations we construct explicit solutions of the nonlinear 
Schrodinger equation with linear and nonlinear periodic potentials. We present ex- 
plicit forms of spatially localized and periodic solutions, and study their properties. 
We put our results in the framework of the exploited perturbation techniques and 
discuss their implications on the properties of associated linear periodic potentials 
and on the possibilities of stabilization of gap solitons using polychromatic lattices. 

Key words: Nonlinear lattices, nonlinear Schrodinger equation, similarity 
transformation, localized modes, Bloch functions. 



1 Introduction 



The significant increase of the interesting properties of the solutions of the 
nonlinear Schrodinger (NLS) equation with a periodic potential during the 
last decade was strongly stimulated by the application of that model in the 
theory of Bose-Einstein condensates (BECs) (see e.g. [Tf2] ). where the model 
is known also as the Gross-Pitaevskii equation. The major advantage of the 
use of optical lattices, which constitute the physical origin of the periodic 
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potential, is the possibility of change of their characteristics in space what 
stimulated studies of managing soliton dynamics by spatially inhomogeneous 
lattices [3J, by low-dimensional lattices [I], by lattices with defects [SHE], etc. 

Independently recently there has been interest on the systematic study of 
nonlinear lattices [8.9. 10J in a more general context, including not only the 
BEC mean-field theory but also electromagnetic wave propagation in layered 
Kerr media. Subsequent studies of the interplay between linear and nonlinear 
lattices [TTIfTS] have revealed more peculiarities of localized modes in such 
structures and unusual stability properties of plane wave solutions. 

Search of the localized modes in all the studies mentioned above was carried 
out within the framework of approximate and numerical methods, as the mod- 
els, in general case, do now allow for exact solutions. However, recently it was 
realized that for a number of NLS models including specific type of periodic 
potentials, exact periodic solutions can be constructed [T3|H] . what can be 
done in a systematic manner, using a kind of inverse engineering and start- 
ing with a given periodic field distribution [lj. We want also to mention the 
Refs. P~51ll6|ll7lfT5| . where different methods to find exact solutions to nonlinear 
Schrodinger equations were studied. 

Even greater progress was achieved in constructing exact solutions, and in 
particular of integrable models, of the NLS equation with time dependent 
coefficients (see e.g. |19ll2Uf21f22ll23f24f25j for the work devoted to equations 
of the NLS-type), and for the stationary GP equation with a potential and 
varying nonlinearity [261127] . However, all mentioned models allowing for ex- 
act solutions look sometimes rather artificial, for example models with exact 
periodic solutions must be infinite or subject to cyclic boundary conditions 
and models with trap potentials require specific laws of time variations of the 
coefficients. Despite those difficulties, they have two important advantages the 
first being that they are still experimentally feasible, and the second being that 
they allow for exact solutions, thus ruling out any ambiguity in interpretation 
of the physical phenomena. 

Returning to the NLS models with spatially periodic coefficients, we observe 
that no exact localized solutions have been reported, so far (the only work to 
our knowledge is the analytical approach for constructing sufficiently narrow 
pulses, recently elaborated in [28]). This is the main goal of this present paper 
to present for the first time such solutions, which can be obtained for specific 
types of linear and nonlinear lattices. Moreover, simple analysis of the pre- 
sented models and their solutions will allow us to make several conclusions 
about properties of the underlying linear lattices and about stabilization of 
localized modes by polychromatic lattices. The results presented in this work, 
represent the extension of the ideas elaborated in earlier publications [261127] 
to the case of periodic nonlinearities. 
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The paper is organized as follows. In Section 2, we develop the theory of 
similarity transformations for our model problem: the nonlinear Schrodinger 
equation with a spatially inhomogeneous nonlinearity. In Section 3, we study 
the stationary localized modes of the inhomogeneous nonlinear Schrodinger 
equation (INLSE). In Section 4, we present a study of the stability of such 
solutions. Finally, in Section 5, we construct explicit periodic solutions of the 
INLSE and study their stability. 



2 Similarity transformations 

We consider the one-dimensional spatially inhomogeneous NLS equation 

iipt = -ipxx + v(x)ip + g{x) \?p\ 2 ip, (1) 

with x G R, v (x) and g(x) being respectively linear and nonlinear periodic 
potentials, whose periods will be required to be equal. More specifically, with- 
out loss of generality we impose the period to be it: i.e. v(x + n) = v(x) and 
g(x + ir) = g(x), what can be always achieved by proper rescaling of variables. 
(It is worth to mention here that in the BEC theory the choice of the period 
tt corresponds to the scaling where the energy is measured in the units of the 
recoil energy). In order to eliminate possible unessential constant energy shifts 
and bring the statement of the problem closer to the standard one we impose 
the requirement for the periodic potential to have a zero mean value, i.e. 

(v(x)) = - r v(x)dx = 0. (2) 

7T JO 

The complex field if)(t, x) will be referred to as a (macroscopic) wave function, 
where again we bear in mind BEC applications. 

In this paper, we focus on stationary solutions of Eq. (flj) , which are of the 
form 

^(t,x)=<P(x)e~^, (3) 

where <j>(x) is a function of x only and \i is a constant below referred to as a 
chemical potential. As it is clear 

-<f) xx + (v(x) - n)<j) + g{x)<f) 3 = 0. (4) 

Now, following the ideas of [26T2T] , we look for a transformation reducing Eq. 
01]) to the solvable stationary NLS equation 

E<& = -$" + G|$| 2 $, (5) 
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where $ = 3>(X), a prime stands for the derivative with respect to X, E and 
G are constants, and X = X(x) is a new spatial variable. To this end, we use 
the ansatz 

0(x)=p(x)$(X), (6) 

where $pf) is a solution of the stationary equation (jSJ) and both p(x) and 
X(x) are functions which must be found from the condition that ip(x) solves 
Eq. (H). 

Substituting (jSJ) into (Hj) we obtain the link: 



= 0, (7) 
as well as expressions for p(a;) and v(x) through p(x) and 



g(x)=G^§-, and i;W = ^+/ 1 -£I I 2 . (8) 
P P 

From the expression for g(x) we find the first constraint of the theory: the 
method is applicable to models with sign definite nonlinearities, and one must 
choose G =sign (g(x)). Next, being interested in solutions existing on the 
whole real axis and excluding singular potentials we have to restrict the anal- 
ysis to functions p{x) in C 2 (lR) which do not acquire zero values and thus 
are sign definite. Moreover since neither (|S]) nor §J§ change as p(x) changes 
the sign (what simply reflects the phase invariance of the model 01])), in what 
follows we restrict the consideration to the case pipe) > 0. 

The solution of (H) is immediate: 



X(x) = / - Tn . (9) 

Jo p {S) 

Here we have taken into account that the constant of the first integration 
with respect to x can be chosen 1, as far as p(x) is still left undefined, and the 
second integration constant conveniently fixes the origin. Thus we have 



X(0) = 0, and lim X(x) = ±oo . (10) 

x— >±oo 



Now the coefficients in Eq. fflh can be rewritten in the form 



(\ Pxx 

E , s G 

x ) = i + V, 9{x) = —. 11 

P P A P 6 
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Thus, ()6]) with (jUJ) gives a solution of Eq. (jlj), which depends on the the 
positive definite function p(x), which must be chosen 7r-periodic. This choice 
is naturally very rich. To be specific we investigate in detail the simplest case 
where 

p(t,x) = 1 + acos(2x), (12) 
with a being a real constant, such that \a\ < 1, and respectively 



4acos(2x) E 
V[X) = ~l + a cos(2x) ~ (1 + a cos(2x)) 4 + ^ ( } 

^ X)= (l + al(2x))6' (M) 
where p is chosen to ensure <^ and reads 

4 2 + 3a 2 

^ 4 -^rrw + E 2(T^r^ (15) 

Some direct generalizations of this model are presented in the Appendix. 

The introduced linear periodic potential v(x) has a number of peculiar prop- 
erties, some of which are: 

(i) It allows for some particular explicit forms of the Bloch function (see e.g. 
expressions (TX5T) and (j2"7j) . below), what is a rather peculiar situation since 
in a general situation explicit forms of solutions of a Hill equation in terms 
of the elementary fucntions are not available. 

(ii) a appears to be the parameter allowing one to control the band spectrum 
(as this is illustrated in each of the panels in Fig. [T]below) and, in particular, 
allowing for smooth transition between continuum spectrum (at a = where 
the periodic potential does not exist and the discrete spectrum (at |a| = 1 
where the periodic potential is represented by a sequence of potential wells 
of infinite depths). In Fig. [1] these two situations are expressed by the fact 
that all gaps are collapsed at a = 0, one the one hand, and one of the bands 
acquire an infinite width at a = 1, on the other hand. 

(iii) By varying the parameter E one can transform the potential from a 
general form with no intervals of coexistence (E < 0) to a form with coexis- 
tence (E > 0) . Hereafter, following [31] under the coexistence phenomenon 
we understand the occurrence when an interval of instability (i.e. a gap) dis- 
appears due to collapsing of two gap edges. This property is also illustrated 
in Fig. [TJ where panels (a) and (b) correspond to the general situation while 
in the panes (c) and (d) one observes intersection of a dashed-doted (red) 
line with the gap edges precisely in the co-existence points. 

It is remarkable that these properties become more clear departing from the 
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solution of the respective nonlinear problem. We consider them in the next 
sections. 



3 Stationary localized modes 

We start with the analysis of the localized modes. To this end we impose 
the boundary conditions lim^-too <f>(x) = and taking into account <^ and 
(1101) we conclude that $(Jf) must satisfy the zero boundary conditions, as 
well: limx^±oo = 0. This is the case where G = — 1 and $(Jf) = 

\J—2E/ cosh(v / — EX) i.e. is the standard stationary NLS soliton. The 

respective solution of Eq. (TjJ reads 



Let us now consider in more detail the obtained solution. First of all we recall 
that the chemical potential of any stationary solution of the NLS equation with 
a periodic potential, which tends to zero or acquires a zero value must belong 
to a forbidden gap of the respective potential and such solutions have a space- 
independent phase, and thus can be chosen real (see e.g. [T1I2U] ). The obtained 
solution (TToT) has chemical potential /x, given by Eq. (FIST) , and hence it must 
belong to a gap of the spectrum of the potential v (x), which is determined by 
the Hill eigenvalue problem: 



Moreover, since the obtained solution (Tl6|) exists for arbitrary negative E, 
including the limit where E — > — oo and consequently fi < v(x) one concludes 
that £ — /j, belongs to the semi-infinite gap (— oo, S[ ) of the band spectrum 
of v(x), i.e. S[ ^ > /j, where we use the notations £^ and £^ for the lower 
and upper edges of the n-th band, i.e. stability region, provided that a = 1 
designates the lowest band. Indeed, assuming the opposite, i.e. assuming that 
at some negative energy the solution belongs to one of the finite gap and 
taking into account the continuous dependence on E one concludes that at 
some negative E, the chemical potential ffl5l) crosses the stability region, what 
contradicts to the existence of the solution with zero boundary conditions. 
Alternatively, the above conclusion follows form the facts that at a = the 
solution (fl6l) is in the semi- infinite gap and that the dependence of /i and v(x) 
on the parameter a G (—1,1) is continuous. 

The described phenomenon is illustrated in Fig. [1] (specifically in the panel 
(a)) where we show the edges of the lowest bands vs the parameter a for four 
typical situations E < (panel (a)); E = (panel (b)); E G (0,E ) with 




(16) 



Vxx + V(x)(p = £(f. 



(17) 
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Fig. 1. Boundaries of the lowest band edges (5 first band edges) of the spectrum of 
Eq. (TTTD as functions of a for (a) E = -5 (b) E = 0, (c) E = 0.1, and (d) E = 1. 

Solid and dashed lines correspond to £n ^ and to £«, . The (red) dashed-dot lines 
represent the chemical potential [i. 

E Q = max a {E m (a)} = 2/5 and E m (a) = 8(1 - « 2 ) 3 /(20 + 15a 2 ) being the 
point of the local minimum of the chemical potential dfi/dE\ E=Eo = for a 
given a, (panel (c)), and E > E (panel (d)). 

Leaving the discussion of the situations depicted in Fig. [1] (c), (d) for Sec. [5j 
now we turn to the limiting cases of the solution (jl~6l) . At E — > — oo it ap- 
proaches the NLS soliton: (j) s (x) ~ v /z 2~E(l + a)/ cosh( v /3 E(l + a)x), what 
reflects the fact that the localization region, determined by 1/ \J\E\, is much 
smaller than the period of the potential. Thus (fl~6]) is an example of an exact 
solution, whose approximate form can be obtained as suggested in [28] (after 
the scaling out the amplitude y/—E). 

An even more interesting situation corresponds to small \E\. When \E\ — > 
the chemical potential approaches the band edge, what is illustrated in Fig. [TJ 
(b), where £ = fi coincides with the edge of the semi-infinite band (in other 
words, here we are dealing with a situation where a formula for dependence of 
the lowest band edge on the parameters of the problem is given explicitly by 
(1151) ). As it is well known (see e.g. [Tf2l?] and references therein) in this case 
the solution is accurately described by the multiple-scale approximation and 
represents an envelope of the Bloch state corresponding to the gap edge. The 
structure of (fl~6l) implies that $(A) is the envelope while 



<p(x) = 1 + acos(2x), (18) 

is the exact Bloch function of the potential v(x, 0), given by (|13p . which cor- 
responds to the lowest edge of the first band. In other words, in the limit 
\E\ — > 0, (jl~6{) is an exact envelope soliton solution bifurcating from the edge 
of the linear spectrum (to the best of authors knowledge (fl~6{) is the first known 
solution of such a type). In Fig. [2] we illustrate the two opposite limits of (fTfil) 
corresponding to an envelope (panel a) and to a narrow (panel b) NLS-type 
soliton. 

Finally we observe that by changing the parameter a, one can scan the semi- 
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Fig. 2. Solutions of Eq. (|4|), which are given by (|16p . for (a) E = —0.01, a = 0.3 
and (b) E = -5, a = 0.1. 

infinite gap, obtaining the localized mode with a priori given detuning towards 
the gap. 



4 Stability of the solutions 



To check the linear stability of the solution (1161) we study the evolution of 
small perturbations of the form ip(x,t) = + f(x,t) + ih(x,t)), where / 

and h are real functions, which leads to the standard linearized Schrodinger 
equation 

/ 




with 



N 



= N 
L_ 

-Li 



(19) 



(20) 



and 



L- = -d xx + v{x) + g{x)(t) 2 {x), 
L + = -d xx + v(x) + 3g(x)(fr 2 (x) 



(21) 
(22) 



For perturbations f,hoc e , we have 



n 2 f = L_L + f. 



(23) 



The operators L_ and L + are self-adjoint. In the following, we will study some 
properties of these operators. Using the fact that satisfies the Eq. PJ, one 
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Fig. 3. [Color online] Stable zones of the solutions (]16p . for different values of the 
parameters a and E (a) -0.5 < E < -0.1 (b) —5<E< -1. 

easily checks that the operator L_ can be rewritten as 

L. = -l 9l (VS, (!■)) . (24) 

As a consequence, / fL_fdx = J \d x \ 2 4> 2 dx > 0, and the operator L_ is 
nonnegative. Thus, the spectrum of the operator L_ is composed of: 

(i) A simple eigenvalue A_ = 0, with the corresponding even eigenf unction, 
which is solution of Eq. (j4]). 

(ii) A strictly positive continuous spectrum [/3,oo). 

With regard to L + operator, it satisfies the following relation L + = L_ + 
2g(x)(p 2 (x). As g(x) is a negative function, the first eigenvalue of L + , A + , 
satisfies A + < 0. So, at least, an eigenvalue of L + is negative. 

Thus, as the L_ operator is nonnegative but the L + operator is nonpositive, 
the composition L_L + is indefinite. Thus, we can not say anything analyt- 
ically of the sign of Q 2 . Moreover, the fact that the potential (JT3l) can take 
both positive and negative values complicates further the analytical study. 
Therefore, we have to resort to numerical methods to solve Eq. f l23l ). 

If some eigenvalue fl 2 is negative, the associated solution is unstable. Other- 
wise, the solution is stable. We can calculate the eigenvalues of the operator 
numerically through a direct discretization of the L_L + operator. In 
Figs. [3] (a) and (b), we show the stable zones of the solutions f|T6|) . for differ- 
ent values of the parameters a and E. In Fig. [3] (a), we have studied the case 
—0.5 < E < —0.1. The zones of stability for the case —5 < E < — 1 are shown 
in Fig. [3]^b). To confirm this result, we have studied numerically the evolution 
of the solutions (|T6|) under finite amplitude perturbations. We have confirmed 
that these obtained solutions are stable in the range given by Figs. [3] (a) and 
(b). 
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5 Exact quasi-periodic solutions 



Let us now turn to (quasi-) periodic solutions of Eq. (0J. To this end, first, 
we have to recall the solutions of the NLS equation §5§ which are expressed 
in terms of the Jacobian elliptic functions, for which we use the standard 
notations [30]. While doing this we bear in mind that different functions $(X) 
which are obtained by the shift of the argument (like, for example, in ( I25a|) 
and f!25dl) below), and thus representing the same solution of the homogeneous 
NLS equation, now will give origin to essentially different solutions of Eq. (J3]), 
because the latter are composed of the two periodic periodic factors steaming 
from the structure of the periodic potential (they are described by the function 
p(x)) and from the NLS solutions mentioned above, (see Eq. (J6])). 

We start with the simplest solutions of Eq. (jlj) for g(x) < (i.e. G = — 1). 
They are 



E 

(f) W (x) = V2ku(l + acos{2x))cn{uX,k), v 2 = - — , (25a) 



• 2) (x) = y/2v(l + ocos(2x))dn(z/X,A;), v 2 = j^-^;, (25b) 

E 



1 -2k 2 ' 
E 

k~2' 



^ 3 \x) = y/2(l - fc>(l + a cobM) k) , ^ 2 = ^, (25c) 

^\x) = y/2(l - k*)ku(l + aco < 2x » 2{lx, % V2= l^2k^ (25d) 

These are real functions, what determines the regions of the parameters where 
they are valid. In particular, taking into account that the elliptic modulus 
k E [0, 1] we have that (x) and (p^ (x) are valid only for E < 0, i.e. these are 
solutions belonging to the semi-infinite gap, while 4>^ l \x) and 4>^\x) belong 
to the gap E < and k > l/y/2 (see Fig. Ufa)) or to the band E > and 
k < l/y/2 (see Figs. (T^c) and (d)). All the solutions bifurcate from the linear 
Bloch state ffl8l) recovered at E — 0. 

The simplest solutions for g(x) > (G — 1) read 



E 

(j) i5) (x) = y/2ku(l + acos(2x))sn(uX,k), v 2 = — , (26a) 

1 + k 

<«) W = V2Ml +a cos(2,))£g|, S-^. (26b) 
As it is clear these solutions are valid for E > 0. 

In the above formulas (I25al) - (126bl) . X(x) is defined in (Tl6|) . and p{x) is given 

by©. 
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Let us now consider the pair of solitons (I26aj) , (126bl) (or alternatively, the pair 
(125IL (|25d])) in the "linear" limit k -> 0. They readily give two eigenstates of 
the potential t> (x) 



<^i = (l + acos(2x))sin(v / ^x) , (27a) 
^ 2 = (l + acos(2x))cos(v A EX) . (27b) 

These are linearly independent solutions. As it is known (see [31 J for the 
details) the coexistence of such solutions occurs only if one of the gaps is 
closed. Thus, in a-dependence of the spectrum of the potential v(x) given by 
(JT5]h f|T5|) at E > there must exits points at which the first lowest gap is 
closes and through which the chemical potential \x given by l fl5l) passes. This 
is exactly what we observed in panels (c) and (d) of Fig. [TJ The existence of 
such points at sufficiently large a steams from the fact that for E > the 
chemical potential infinitely grows with a while the minima of the periodic 
potential tend to — oo. 

Other exact eigenstates of the potential v(x) given by (fl3|) can be obtained by 
considering the linear limit of the solutions solutions f!25cj) and (125dj) which 



corresponds to k — > 1. In this way we obtain a pair of two "unstable" solutions 



^ 3 = (l + acos(2x))cosh(v /3 EX) , (28a) 
</? 4 = (l + acos(2x))sinh(v /z ^x) • (28b) 

As it is clear the combinations ip 3 ± ip^ supplied by the proper constant factor 
provide the asymptotics of the solitary wave 4> s (x) given by (1TB]) at x — > =Foo. 

Taking into account that g(x) is sign definite, simple stability analysis of the 
presented solutions can be performed following |14] : since <f>W(x) > and 
^ 3 \x) > one verifies that they are linearly unstable, while the stability of 
the solutions (f)^(x), (f)^'(x), <p^\x) and (j)^(x) is left undetermined. 

Finally, we consider the limit k — > 1/V2 of the solutions (I25a[) and ( I25dl) . 
which corresponds to the limit of large nonlinearity. For this case, the potential 
v(x) can be viewed as a small and smooth perturbation of the NLS equation, 
whose nonlinearity is also a slow function of the spatial variable. Thus the 
mentioned solutions can be viewed as the periodic NLS solutions modulated 
by the "envelope" <po(x)- 
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6 Physical applications and concluding remarks. 



Let us now turn to the physical applications of the obtained results. This 
problem arises naturally since the periodic structures we have used, Eqs. ( TTBl 
and ffl~4"j) . being even expressed in elementary functions, are still not feasible 
for the most experimental settings, where only one or a few laser beams are 
used (if one bears in mind optical lattices for BEC applications). One thus can 
pose the question as follows: do the obtained solutions represent satisfactory 
approximations to some realistic localized modes (like for example the ones 
found numerically in [11]) in the models where the periodic coefficients are 
represented by a few first Fourier harmonics of the potential v(x,E)7 The 
present section aims to answer this question. 

To this end we Fourier expand the functions v (x) and g(x), and introduce 
the truncated potentials Vn(x) and Qn^x) generated by superpositions of iV 
harmonics: 



N N 

Vn{x) = E ) cos(2ra), g N (x) = 9n(a) cos(2nx), (29) 

n=l ra=0 

where, for N = 2 



/ t I \ „1 - vl - a 2 «(4 + a 2 ) ,„„ s 



(1 - Vl - a 2 ) 2 5a 2 
^ E ) = L - E {1 - a 2y,2 > ( 30b ) 



and 



. . 140a 2 + 8 + 15a 4 /oi . 

9o(u) = -- (1 _ a 2)ii/2 » ( 31a ) 

3a 8 + 12a 2 + a 4 

9M) = ^ (1 _ a2) ii/2 » ( 31b ) 
21a 2 2 + a 2 

92( a ) = ~ 4 (1 _ a 2)ii/2 - ( 31c ) 



Next we use the Eqs. ( 129|) . (|30aj) . (I30bl) and fl31al) . (I31b|) . (I31cj) to approximate 



the functions f (x) and ^(x), given by Eqs. (|T4|) and (|T3|) . respectively. Thus, 
for example, for a = —0.1 and i£ = —0.5, which are values where the solutions 
of Eq. (@J are stable, (see Fig. G^a)), we obtain 
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Fig. 4. Evolution of the solutions CEj), with x G [-40, 40], t € [0, 1500], for £ = -0.5 
and a = —0.1 (a) using one harmonic (b) two harmonics. For the case (a), the 
solution is oscillatory. For the case (b) the solution is stable. 



Finally, we have simulated numerically the dynamics of the wave packet with 
the initial profile given by the solution ( jl~6l) . and described by the evolution 
equation (CQ) with the potentials vn and qn with N = 1, 2 instead of v and 
g, respectively. Some typical results are shown in Fig. HI One observes that 
while the explicit analytical solution ({TBI) is not a satisfactory approximation 
for the harmonic potentials: in Fig. H](a) one observes oscillatory behavior of 
the mode. Using a two-harmonic approximation for the potential results in a 
very stable behavior of the solution. 

It turns out, that for some specific domains of the parameters the obtained 
solutions represent fairly good approximations even for the harmonic poten- 
tials. An example of such a situation is shown in Fig. [5] where we have choses 
E = —0.1 and a = —0.1: no visible difference exist in the dynamical regime 
for the monochromatic lattices (Fig. [5](a)) and the lattices in a form of super- 
position of the two harmonics (Fig. 0(b)). 

Both Figs. H] and [5] were obtained by direct numerical simulations of Eq. (CQ), 
using as initial data in the evolution in time the solution (fTBl) . 

To conclude, in this paper, by using similarity transformations we have con- 
structed explicit localized solutions of the nonlinear Schrodinger equation with 
linear periodic potential and spatially periodic nonlinearities. We have studied 
such solutions and their linear stability. We also have calculated periodic so- 
lutions of the inhomogeneous nonlinear Schrodinger equation and have shown 
that they reveal some properties of the underline linear lattices, providing 



g = -1.1099 - 0.6436 cos(2x) - 0.1115 cos(4x) 



(32a) 
(32b) 
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Fig. 5. Evolution of the solutions ([IB]) , with x G [-40, 40], t £ [0, 1500], for £ = -0.1 
and a = —0.1 where the linear and nonlinear lattices are approximated by using 
only one harmonic (a) and two harmonics (b). For both cases, the solution is stable. 



exact analytical expressions for the respective Bloch waves. 
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A Examples of models allowing exact solutions 

In this appendix, we present two more general models that we used in Eqs. 
(TTBl and (|T4|) . The first model corresponds to the choice 



p(x) = (l + acos(2x)) p 



(A.l) 



with peR. So, it is clear that 



G 



(A.2) 



(1 + acos(2a;)) 6 P' 
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and the external potential becomes 



v(x) = —Aap (1 + a cos(2x)) a(l — p) + cos(2x) + pa cos 2 (2x) 

E 



+ fi. (A.3) 



(1 + acos(2x)) 4 P 

It is clear that when p = 1, we recover the expressions ( jl~3j) and (JHJ). 

The second model corresponds to 

g(x)= V (dn(ti,k)) p , (A.4) 

with 77 a constant and dn the Jacobi elliptic function. We can calculate the 
external potential v(x) and obtain 



«(*) = |[(l-f) (dn(^)) 2 ]+^ [(^ 2 -l) (l + f) 





fc 2 y 


ft- 


6 J. 



67 (dn(£,fc)) 2 
-E(|) 2 / 3 (dnK,fc)) 2P/3 + ^ (A.5) 

One can obtain the nonlinearity g(x) and the potential v(x), according to the 
value of p. 
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